Strongly Correlated Electrons in the [Ni(hmp)(ROH)X] 4 Single Molecule Magnet: A DFT+U 

Study 
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The single-molecule magnet [Ni(hmp)(MeOH)Cl] 4 is studied using both density functional theory (DFT) 
and the DFT+U method, and the results are compared. By incorporating a Hubbard-U like term for both the 
nickel and oxygen atoms, the experimentally determined ground state is successfully obtained, and the exchange 
coupling constants derived from the DFT+U calculation agree with experiment very well. The results show that 
the nickel 3d and oxygen 2p electrons in this molecule are strongly correlated, and thus the inclusion of on-site 
Coulomb energies is crucial to obtaining the correct results. 

PACS numbers: 75.50.Xx, 75.30.Gw, 71.I5.Mb, 71.10.Fd, 73.20.At 
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Single-molecule magnets (SMMs) have drawn much atten- 
tion since their discovery in 1991 til S 13J. SMM crystals 
contain ordered arrays of molecular nanomagnets, each pos- 
sessing a large spin ground state (S = 10 for Mni2-Ac) and a 
significant uniaxial magneto-anisotropy (DS 2 , with D < 0). 
These two ingredients give rise to a magnetic spectrum for 
an isolated molecule in which the lowest lying levels corre- 
spond to the 'spin-up' and 'spin-down' states (m s = ±5*), 
separated by an energy barrier of order DS 2 . This barrier re- 
sults in magnetic bistability and hysteresis at low temperatures 
(ksT << DS 2 ). In contrast to bulk ferromagnets, however, 
this hysteresis is intrinsic to the individual molecules. There 
has, therefore, been much interest in the potential implemen- 
tation of SMMs as the elementary memory units in both clas- 
sical and quantum computers. 

For most transition metal complexes (including Mni2-Ac), 
the intramolecular superexchange between the constituent 
ions is predominantly antiferromagnetic (AFM). Neverthe- 
less, due to spin frustration effects, uncompensated moments 
of many tens of [Ib are often realized. However, the ability to 
engineer pure ferromagnetic superexchange within a molecule 
is highly desirable J3l, because this ultimately removes one 
of the many challenges in designing new and better SMMs. 
A relatively new series which has attracted recent interest is 
[Ni(hmp)(ROH)Cl] 4 fi @, 0, % % (hereon denoted Ni 4 ), 
where hmp is the anion of 2-hydroxymethylpyridine, and R is 
an alkyl substituent such as methyl, ethyl, etc. Several exper- 
iments, including EPR studies |6] and magnetic susceptibility 
measurements ||5|] clearly show that the ground state of Ni 4 is 
ferromagnetic with total spin S = 4. In this letter, we present 
the results of detailed density functional theory (DFT) calcu- 
lations (including on-site Coulomb energies) which provide 
crucial insights into the origin of this ferromagnetic state. 

While DFT ifioll has successfully explained the proper- 
ties of a variety of SMMs, including Mni2, M114, C04, Fe4 
O El El Hill, and even some other nickel based SMMs 
01611 . it has so far failed miserably for N14. Not only were 
the early theoretical attempts unable to reproduce the correct 
ground state, but the resulting coupling constants were also 
found to be antiferromagnetic, and orders of magnitude higher 




FIG. 1: Optimized structure (a) and exchange-coupling scheme (b) 
of [Ni(hmp)(MeOH)Cl] 4 . Magenta atoms are Ni; red atoms O(l); 
yellow atoms CI; orange atoms other oxygens (i.e. 0(2)); blue atoms 
N; and gray atoms C. The Ni and O(l) atoms form a slightly distorted 
cube. 



than the experimental values fl7ll . It was also found in the cal- 
culation that the spin density is not quite localized around the 
nickel atoms, as expected. Thus, it has been suggested that the 
discrepancy between theory and experiments might arise due 
to the small "spin density leakage" in this system, resulting in 
spin derealization. 

There is another possibility. Due to the localized nature of 
3d electrons, transition metal dioxides, including nickel ox- 
ides, are known to be strongly correlated materials. The func- 
tioning core of [Ni(hmp)(ROH)Cl] 4 , on the other hand, is a 
cubic tetra-nickel oxide (N14O4), which is structurally very 
close to the nickel oxide complex. Therefore, it is more prob- 
able that the lack of strong correlation in DFT is responsible 
for this failure, and the "spin density leakage" is just an arti- 
fact. To justify this speculation, we calculated the electronic 
structure of [Ni(hmp)(MeOH)Cl] 4 using both the DFT and 
DFT+U methods. The latter was introduced by V. I. Anisi- 
mov et al. 01811 and simplified by M. Cococcioni et al. IU9I1 . 

All the reported calculations were done using the PWSCF 
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package |20ll . which utilizes PBE exchange-correlation func- 
tionals [21], ultrasoft pseudopotentials 12211 . and aplane-wave 
basis-set. We respectively chose energy cut-offs for the wave 
functions and charge densities to be 40 Ry and 400 Ry to en- 
sure total energy convergence. The structure of the molecule 
was optimized with a fixed total spin 5 = 4 until the force 
on each atom was smaller than 0.01 eV/A. The relaxed struc- 
ture is in good agreement with experimental results. The same 
structure was then used for some of the AFM states (S = 0) 
and S = 2, as well as for all of the DFT+U calculations. Due 
to symmetry restrictions, we only simulated the 5 = 4, 5 = 2 
and AFM (5 = 0) states. For the DFT+U calculations, a self- 
consistent Hubbard-U method [19] has been incorporated to 
determine the U value for Ni, which turns out to be 6.20 eV 
for this system. For oxygen, we took the well-established U 
value of 5.90 eV H. For the DOS and projected DOS, we 
used 0.1 eV gaussian smearing to smooth the results. 

The DFT calculations confirmed Park et al.'s findings. The 
optimized structure is shown in Figure [T(a)| In order to better 
display the geometry, we hide all hydrogen atoms. The four 
nickel atoms and four oxygen atoms on the hmp group (we 
call them O(l) from now on) define a slightly distorted cube. 
The AFM state (S = 0) turns out to be the ground state, which 
is 14.8 meV lower than the 5 = 2 state and 35.1 meV lower 
than the 5 = 4 state (table |U column DFT). Using Lowdin's 
charge analysis, one can see that about a 1.48-1.50 fiB mag- 
netic moment is found on each nickel atom, and each of the 
four O(l) atoms contributes a 0.1-0.26 fiB magnetic moment 
(table [Til column DFT). The Heisenberg Hamiltonian in gen- 
eral can be written as: 



Hex — ^ ' Jij^i ' Sj 



(1) 



Considering the molecule's 54 symmetry (Figure |l(b)[ ), the 
six JijS reduced to just two values, i.e.: 



Hex — Ji(Si ■ S2 



S2 ■ S3 + S3 • S4 + S4 ■ Si) 

+J 2 (S! • S 3 + S 2 • S 4 ). (2) 



We then determined these exchange coupling constants by fit- 
ting the expression of Eq. (2) to the obtained energies of 
the 5 = 4, 2 and spin states, giving J\ = 3.54 meV and 
J2 = 2.12 meV. The positive sign for both numbers indicates 
AFM coupling. From the total electronic density of states 
(DOS), as well as the projected density of states (PDOS) onto 
Ni (Figure|2]i, it is clear that, in the DFT calculation, the nickel 
atom 3d orbitals dominate both the highest occupied (HOMO) 
and lowest unoccupied molecular orbitals (LUMO), so that 
the HOMO-LUMO gap is of the d-d type. 

The DFT+U calculations were performed in two stages. In 
the first stage (we call this DFT + U d ), we turned on the 
Hubbard-U parameter for the nickel 3d orbitals only, like 
most current DFT+U calculations; in the second stage (we 
call this DFT + U p+d ), we turned on the U parameter for 
both the nickel 3d and O(l) 2p orbitals. In fact, it is known 
that Coulomb interactions between oxygen 2p electrons are 
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FIG. 2: Total DOS and PDOS on the 2p and 3d orbitals of 
[Ni(hmp)(MeOH)Cl] 4 for the S = 4 state using DFT. Red lines 
represent a-spin and blue lines represent /3-spin. The green line is 
the Fermi level, Ef- The DOS were drawn to the same scale for 
comparison purposes. 





DFT 


DFT + U d 


DFT + U p+d 


AFM (S=0) 


0.0000 


0.00000 


0.000000 


S=2 


0.0011 


0.00012 


-0.000069 


S=4 


0.0026 


0.00019 


-0.000368 



TABLE I: Total energies in Rydbergs. All numbers are relative to the 
AFM state (5 = 0) 



comparable to those between d electrons 11241 12511 . and should 
hence be taken into consideration as well. However, since 
oxygen usually bares a fully occupied p-shell, this correlation 
effect is often thought to be neglegible. Therefore, in most 
cases, DFT + U d can already yield a satisfactory description 
of the ground-state without oxygen 2p-electron corrections. 
Nevertheless, DFT+U has to be taken into consideration ex- 
plicitly here for both the 3d and oxygen 2p electrons in order 
to obtain the correct ground state for this molecule. 

The DFT+U energies are shown in table Q] By turning on 
DFT+U for the nickel atoms only, the energy differences be- 
tween different spin states were greatly reduced, hence giv- 
ing much smaller numbers for the exchange coupling con- 
stants. However, the ground state here is still AFM (5 = 0), 
and the energies for the 5 = 4 and 5 = 2 states relative 
to the AFM state are 2.61 meV and 1.60 meV, respectively. 
But once we take into consideration the strong coulomb in- 
teractions for both the Ni and O(l) atoms, the order is re- 
versed, yielding correctly a 5 = 4 ground state. A spin- 
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FIG. 3: Total DOS and PDOS on the 2p and 3d orbitals of 
[Ni(hmp)(MeOH)Cl] 4 for the 5" = 4 state using DFT + U d . Red 
lines represent a-spin and blue lines represent /3-spin. The green 
line is the Fermi level, Ef- DOS were drawn to the same scale for 
comparison purposes. 
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FIG. 4: Total DOS and PDOS on the 2p and 3d orbitals 
of [Ni(hmp)(MeOH)Cl] 4 for the S = 4 ground-state using 
DFT + U p+d . Red lines represent a-spin and blue lines represent 
/3-spin. The green line is the Fermi level, Ef- DOS were drawn to 
the same scale for comparison purposes. 







DFT 




DFT + U 


d 


DFT + U p+d 




AFM 


S=2 


S=4 


AFM 


S=2 


S=4 


AFM 


S=2 S=4 


Ni 


1.49 


1.49 


1.50 


1.68 


1.68 


1.68 


1.74 


1.74 1.74 


0(1) 


0.11 


0.13 


0.26 


0.06 


0.08 


0.16 


0.04 


0.06 0.11 


CI 


0.09 


0.09 


0.09 


0.06 


0.06 


0.06 


0.06 


0.06 0.06 


N 


0.08 


0.08 


0.08 


0.06 


0.05 


0.05 


0.05 


0.05 0.05 


0(2) 


0.06 


0.06 


0.06 


0.03 


0.03 


0.03 


0.03 


0.03 0.03 



TABLE II: Magnetic moments (in /is) captured by Ni, O(l), CI, N 
and 0(2) atoms. AFM indicates the antiferromagnetic state (S = 0). 
All numbers are averaged over the same species. 



unrestricted calculation also confirmed this discovery. The 
S — 2 state is now 0.94 meV lower than the AFM state, and 
the S = 4 ground state is 5.00 meV lower. Using these val- 
ues, we obtained ferromagnetic exchange-coupling constants 
for the DFT + U p+d calculation from a fit to equation 2, i.e. 
Ji = —0.50 meV and J2 = —0.68 meV|2] These results 
match experiment reasonably well (—0.68 and —2.28 meV) 
& 




Up (eV) 

FIG. 5: Variation of Ji and J2 with respect to different U p for O(l). 



To better understand the contribution of the Hubbard-U like 
term, we first performed Lowdin's charge analysis to calculate 
the magnetic moments captured by the Ni, O(l), CI and N and 
0(2) atoms (table HJi. The results show that the spin density 
is more localized in the DFT + U p+d calculations (1.74 /xg) 
than in normal DFT (around 1.50 /x^), and that the magnetic 
moments found on the O(l) atoms are greatly reduced. This 
is because the strong on-site Coulomb interaction prevents the 
hybridization between the nickel 3d and oxygen 2p orbitals, 
thus preventing the unphysical "spin-leakage". The DFT cal- 
culations favor the AFM ground state because the lack of on- 
site energy tends to couple electrons with opposite spin pro- 
jections, and thus lead to the incorrect ground state. The lo- 
cal magnetic moment on an individual atom is a measurable 
quantity using NMR, and thus can be used to validate these 
theoretical predictions. 

Total DOS and PDOS in the DFT + U d method (Figure 
[SJand DFT + U p+d method (Figure |4|i were also calculated 
H27I1 . In contrast to the DFT results, the dominant contribution 
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to the HOMO and LUMO states in both DFT+U calculations 
is now from the 2p (CI and O) and 3d (Ni) orbitals, respec- 
tively (FigureHJ. We do not have direct experimental results to 
compare with this feature, however, for nickel oxide it is well 
known that DFT gives incorrect PDOS contributions. Early 
experiments and calculations 1 28 , 29i 30l 31, 32] show that 
instead of a d-d gap given by DFT, nickel oxides actually have 
a p-d gap. The largest spin density contribution, of course, is 
still from the Ni 3d electrons. For the S = 4 ground state, 
the DFT + U p+d calculation yields a LUMO-HOMO gap of 
2.95 eV, which is from majority spin to minority spin; in the 
DFT + U d and DFT calculations, these numbers are 2.56 eV 
and 1 .09 eV, respectively. Since DFT has been known to un- 
derestimate energy gaps and excitation states, DFT+U calcu- 
lations have proven necessary in order to obtain agreement 
with experiments such as resonant inelastic X-ray scattering 
(RIXS) and XPS [133]. The present study also likely calls into 
question recent reports of HOMO-LUMO gaps of fractions of 
an eV in the Fe 8 SMM Q- 

Finally, we have also repeated DFT + U d and 
DFT + U p+d calculations with U d for nickel ranging 
from 4.58 to 6.20 eV and U p for oxygen ranging from 1.0 to 
8.0 eV. The energetic order of spin states remains the same, 
and the magnitudes of the energy differences between spin 
states are more-or-less insensitive to variations of U p from 
3.0 to 8.0 eV, as seen in Figure [5] This clearly demonstrates 
the reliability and robustness of our results. 

In conclusion, we have performed DFT and DFT+U cal- 
culations for [Ni(hmp)(McOH)Cl] 4 . Because of the strong 
correlation effects in this system, the DFT calculation fails 
due to the fact that the lack of on-site energy unphysically 
encourages the hybridization of orbitals, leading to AFM cou- 
pling. The inclusion of a Hubbard-U like term for both the 
Ni 3d and O(l) 2p electrons greatly enhances the localization 
for both states, and is essential in order to obtain the correct 
ferromagnetic ground state and exchange-coupling constants. 
After taking both corrections into consideration, these prop- 
erties were successfully reproduced by the calculations. We 
then analyzed the DOS and projected DOS of the system, and 
the calculation predicts that the optical transition from HOMO 
to LUMO is p-d like, and the gap is 2.95 eV. 
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